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TECHNICAL PAPER 


PLATE AND BUTT-WELD STRESSES BEYOND ELASTIC LIMIT, 
MATERIAL AND STRUCTURAL MODELING 

I. INTRODUCTION 


The solid rocket booster (SRB) aft skirt is a 22I9-T87 aluminum, high performance structure that 
supports the entire shuttle vehicle through stacking operations into liftoff. Two flight copies were tested 
to revised operational loads, and both failed at the same weld location at considerably less load than the 
estimated ultimate. Monitoring gauges also indicated the failure strain to be much less than predicted. 
This modeling and test prediction breakdown provoked an effort to understand the related physics ot 
stresses beyond the elastic limit. 

What makes the inquiry of inelastic behavior so compelling is that all high performance structures 
must embody this phenomenon in the ultimate safety factor analyses. Forty percent of the alt skirt 
strength is sustained by material properties beyond the elastic limit, a region not too well understood. 

To gain basic insights on inelastic structural behavior, an analytical modeling approach was 
developed which incorporates nonlinear material properties with familiar linear strength ol materials 
method. Material properties through the inelastic range were modeled from uniaxial test data and were 
applied to analysis of rectangular bars and plates subjected to tension, bending, and combined loadings. 
Associated stress and strain equations were derived in sufficient detail for the collective understanding 
and benefit of all interfacing disciplines, especially design, stress, materials, and test. 

The analytical method is applicable to most strength of materials problems and may be extended 
to some classical elasticity models. It should be especially applicable to combustion device structures and 
aeroheated zones, where high temperatures decrease the elastic limit while extending the inelastic range 
of materials. 

Metallurgical discontinuity stresses resulting from dissimilarity of weld and parent material 
properties beyond the elastic region are included. A previous paper [ I ] explored this phenomenon and in 
this report is made consistent with the succeeding material and analytical modeling technique. 

Of 24 structural failures in Marshall Space Flight Center (MSFC) programs, the average ultimate 
test loading fell 15 percent below the predicted safety factor [21. It is hoped that this paper will help to 
better understand the inelastic failure of elementary structures through the analytical method, will 
provide a quick backup analysis for comparing with new and existing computational methods, will 
identify practical limitations, and will support safe design analyses and test data evaluations. 



II. MATERIALS MODELING 


Curve OACE of figure 1 typifies the stress-strain data obtainable from a uniaxial tensile test of a 
polycrystalline material. Points A and E are the elastic limit and the ultimate stress properties, respec- 
tively. Line OA is the elastic portion of the curve which depends on the resilience between molecules. 
Plastic flow is the permanent deformation caused from displacement of atoms, or molecules, to new crys- 
tal lattice sites. When the applied stress exceeds the elastic limit, elastic straining decreases while plastic 
flow rate increases. The ratio of elastic strain and plastic flow defines the slope, and their change rate 
characterizes the nonlinear property of curve ACE as it approaches full plastic state. 

The area under the stress-strain curve OACE represents the total strain energy per unit volume of 
the material. When the material is stressed inelastically to point A' and relaxed, it relaxes elastically 
(same slope as OA) to a new origin point O’ for total plastic strain of O-O' . The material takes on a 
permanent set, and the strain energy represented by the area O AC A' O' is not recoverable. Consequently, 
the new stress-strain property is defined by O'A'E , and the elastic limit advances to point A' . The plastic 
slip along the crystal boundary is directional, which transforms the isotropic material to anisotropic in the 
next load cycle. 

Current structural design practice is based on linear material properties in which the stress-strain 
relationship to ultimate is defined by the slope of line OA extended to point B. Point B is equal to the 
ultimate stress at point E. This assumption is rather conservative for local discontinuity stresses because 
it embraces a small fraction of the total strain energy available for the ultimate safety factor. When opera- 
tional loads on inventory structures are projected to increase beyond original design loads, it is prudent to 
calculate a new safety factor using the total strain energy, in which the increased energy may be sufficient 
to establish that actual margins are still within specifications. 

The success of a nonlinear determination depends on the accuracy of the finite element method 
(FEM) code and structural modeling, and on the understanding of results. Sometimes, stress results of 
local areas from large commercial codes are so intricate and the mechanics so obscured that a hand analy- 
sis is necessary to interpret the physics. Often, new codes and options are difficult to verify, and they 
succumb to cookbook application. When this happens, a backup hand analysis of basic structural 
elements for comparison and perception becomes essential, and its applicability depends on the sim- 
plicity and compatibility of inelastic material expressions. 



Strain 

Figure I. Design data from uniaxial test. 


2 


A. Uniaxial Stress-Strain 


Uniaxial tension tests are the simplest type for obtaining mechanical properties of structural 
materials. Stress-strain material properties of commonly milled bars and plates are well defined and 
widely published [3], but properties of specially processed materials, such as welds, castings, and forg- 
ing, are usually developed by the user. Such was the case of the SRB skirt butt-weld which is a good 
basis for this material modeling discussion. Properties defined and qualifications expressed are not 
standard, but may provide guidelines for designing material model tests applicable to analytical 
approaches. 

Multipass butt-weld properties vary significantly with design geometry which influences the heat 
intensity and distribution across the weld and into adjacent heat affected zones (HAZ). Obviously, the 
stress-strain relationship across the weld width varies uniquely for each design and manufacturing 
process and must be obtained from sample test specimen. If properties of the HAZ are required and the 
width of the zone is less than the gauge length, similar properties in the adjacent overlapped zones must 
be determined and separated from the gauge data through the rule of mixtures. All this points to the dif- 
ficulty of not only obtaining weld properties, but the necessity of discussing the data with the intended 
user. 


The normal stress-strain test data for the SRB skirt analysis, presented in figure 2, was developed 
from 24 specimens [4], using O.I25-in electrical strain gauges mounted at midwidth on weld surfaces. 
The data is the statistical average of the sample size and is preferred for analytical methods that are to be 
used for tracking test instrumentation. Deviation of the data about the average is used to define the 
A-basis ultimate properties required for design analysis and test margin predictions. 



Since stress is related to strain through material properties, and as stress increases beyond the 
elastic limit, these properties change nonlinearly with stress. The simplest normal stress-strain model 
representing these changes is a two-parameter curve-fit. 


a 


ae 


b 


(I) 
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where a is the strength coefficient and b is the strain hardening exponent. When b = I , the expression 
reduces to Hooke’s law of elasticity, and coefficients is identically Young's modulus, E. When b = 0, 
pure plastic flow occurs at constant stress, and the material is said to be perfectly plastic (appendix C). 

The two parameter model is particularly compatible with the analytical method because material 
properties expressed as coefficients and exponents of stress or strain are easily combined with stress and 
strain variables of strength formulas through their coefficients and exponents. 

A conservative material representation should curve-fit the origin, elastic limit, and ultimate 
stress parameters without exceeding the area under the test data curve, as shown in figure 2. Data 
required to determine the two parameters of equation ( 1 ) are specified by properties at the two extreme 
conditions of the inelastic curve: Young's modulus, E ; elastic limit stress, ay/ ; ultimate stress, oy,-: and 
ultimate strain, e t /. The elastic limit stress, though difficult to isolate, is a more meaningful property in 
the inelastic process than the yield point. If maximum stress of very ductile materials precedes fracture, 
strains beyond maximum stress should be considered unstable, and the value at inception of maximum 
stress is assumed ultimate stress. 

A better midrange curve-fit would be possible by using the true strain rather than the suggested 
engineering strain. However, the slight gain in accuracy would severely compromise the potential appli- 
cation of the analytical process. Besides, the error decreases as the stress approaches failure. Since proto- 
type test predictions require accurate midrange properties, a second set of materials models should be 
derived between the elastic limit to 20 percent beyond. At the same time, effects of maximum expected 
error should be weighed against effects of the coefficient-of-variation of material data which has been 
reported to exceed 20 percent |4] . 

The hardening exponent of equation (1) may be calculated from 


b = log(a/;/CT/r / )/log(e £ y*£/a / r / ) 


( 2 ) 


and the coefficient from 


a = (Tf;,J((T, : ,/E) h 


(3) 


Subsequent material models derived from simple uniaxial tension data are expressed in stress or 
strain variables and cast in similar two-parameter form. This is a necessary condition for simplifying the 
analytical process. 

An A-basis strength of material is defined as having a 99-percent reliability taken at a 95-percent 
confidence level. This tolerance limit is used for both elastic limit and ultimate strength. The A-basis 
strength is calculated from referenced weld test data having average strength of a^ = 48.6 ksi, standard 
deviation of s = 1.6 ksi, and a factor of k = 3.18 for a 24 sample size. Applying these values to the 
statistical criteria, 
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(J(J = (T — ks , 


( 4 ) 


results in a design ultimate stress of <t (! = 43.5 ksi. The A-basis elastic limit stress is o>/ = 20.5 ksi. As 
the number of test specimen increases, the A:- factor decreases to the limit of 2.33. which has trades 
implications as to the number of specimen required versus performance of the structure. 

It must be pointed out that the material model specified above is compromised in three respects. 
First, equation (1) assumes a nonlinear stress-strain relationship throughout. This analytical method is 
definitely not applicable below the elastic limit, except as required to derive inelastic bending models. A 
more accurate relationship would have split the curve-fit into two parts; a standard linear part having b = 

I and a = E up to the elastic limit of 26.3 ksi followed by the nonlinear part. However, incorporating 
this discontinuity in bending stress derivations would be extraneously complicated and will be shown to 
be insignificant beyond the elastic limit for which it is intended. 

Secondly, the stress-strain curve of figure 2 assumes a symmetric material. There is reason to 
believe that most crystalline materials have different tension and compression elastic limits which may be 
demonstrated in a pure bending experiment. However, this is an order of refinement which is not neces- 
sary with current general developments. 

Finally, it should be recognized that weld properties are not only manufacturing process 
dependent, but are also dependent on the mass distribution of the adjoining material along the weld 
length. If the weld is to join a highly sculptured forging, such as the skirt post, the heat sink will vary 
accordingly, as will the weld and the HAZ properties. Then the materials representation and number of 
models to be developed must be decided jointly by manufacturing and stress considerations. Does the 
weld thickness, therefore weld-heat, taper along the length? Does the maximum stress and weakest 
expected weld occur at independent places along the weld length, etc.? Similar questions may be asked 
of forgings and castings properties. 


B. Poisson’s Ratio 

All polycrystalline materials possess two elastic constants from which all other constants may be 
derived. Young’s modulus and Poisson’s ratio may be obtained simultaneously from a simple uniaxial 
test specimen. Poisson’s ratio is defined as the ratio of unit lateral contraction and unit axial extension, |x 
= —e } Je x — — s-/e v , which suggests a volumetric relationship. When a material in uniaxial test is 
extended, the space between molecules in the crystal lattice are displaced linearly causing a minute 
volume increase. The new volume is related to the specimen dimensions and strains by 


V = (L + Le x )(t + te v )(B + Be.) . 


Substituting Poisson's ratio expressions for lateral strains and ignoring high order terms, the volume 
change relates to Poisson’s ratio by 
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In the elastic range, Poisson's ratio is a constant having a value between 0.28 to 0.33 for most 
common materials, and 0.3 is used in most general analyses. When a specimen is extended beyond the 
elastic limit, molecules that slip into new lattice sites cause no change in volume. Therefore, a totally 
plastic Poisson’s ratio engages in no volume change and should not exceed fjq, = 0.5. Consequently, as 
the axial strain increases beyond the elastic limit, the volume change rate decreases and Poisson’s ratio 
approaches its maximum value. 

Because the same molecular displacement process which occurs in the normal stress-strain 
relationship also is volume related, it would seem reasonable to expect Poisson’s ratio to be expressed 
similarly to equation (1), 


P = ge* ■ 


(5) 


Current practice is to express the Poisson’s ratio nonlinearly with the secant modulus defined by 
normal stress-strain data. Using data developed in figure 2 and equation (1), Poisson's ratio based on 
secant modulus is 

I lb 

P = P/> [ P/4 Pf 1 y |^-J (6) 

where subscripts p and e refer to plastic and elastic properties, respectively, and E is Young’s modulus. 
This expression is adjusted to match the elastic Poisson’s ratio with the elastic limit stress v FL through the 
variable given by. 



(7) 


Equation (6) is incompatible with the subsequent analytical method for reasons given in the above 
stress-strain expression. However, the secant model is still prerequisite to calculate the Poisson’s ratios at 
the elastic limit stress p, FL and at the ultimate stress p.^. Using stresses from equation (1) and Poisson’s 
ratio from equation (6), parameters of equation (5) are calculated similarly to the tension model: 


/ = logtixy/pL/r^j/logla^/a^] 1 ^ , 


( 8 ) 


8 — P ■Ed(®Ed a 'f b ■ 


( 9 ) 



Figure 3 presents the two constructed Poisson’s ratio models discussed, based on uniaxial test 
data of figure 2. Though common end conditions were imposed on both models, a deviation of up to 10 
percent is noted at midrange between the two plots. This deviation may affect monitoring and evaluating 
strain gauges on prototype structures. How much uncertainty is tolerated may be determined from sensi- 
tivity analyses of specific structural environments. Poisson's ratio expressed by equation (5) is more 
linearly related to stress and provides another option to correlate with actual test data. 



Strain, in/in 


Figure 3. Constructed Poisson’s ratio versus strain. 


C. Failure Criteria 

The total strain energy produced in multistress fields consists of energy causing a volume change 
and energy causing distortion. Basis for the failure theory is that distortion energy is limited and that 
hydrostatic strain does not cause failure. Therefore, hydrostatic strain is deleted from the total strain 
energy [5]. The distortion energy theory of failure hypothesizes that a material will fail when the energy 
of distortion of combined stresses achieves the same energy for failure in simple tension. 

Using a two-dimensional stress field, the failure criterion is expressed by 


°t/ ^ crr + CTi - CT l°2 . 


( 10 ) 


where subscripts 1 and 2 refer to principal stresses. Roots of the characteristic equation formed from the 
two-dimensional stress tensor. 


(cr tr -X) t,, 
T„ (cr-, - X) 


= 0 , 


define the inplane principal stresses, 
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__ + cr ; 1 


CT 1.2 = V (a r + o\) 2 — 4 (ct v ct- — t v :) 


Substituting these roots into equation (10) gives the biaxial stress failure condition, 


c Ju 5 Va; + ct- 2 - ct ( ct ; + 3 t a ? . 


(ID 


Shear stress is noted by t, and u v is the maximum stress obtained from a uniaxial tension test. Note that 
when biaxial normal stresses have unlike polarity, the induced distortion energy to failure increases. 
When stresses have the same polarity, distortion energy decreases and is a capital benefit in pressure 
vessel designs. 

The three-dimensional failure criteria is 


&U 


+ 02 + 0"3 O' I <*2 — <T|03 — OSO^ , 


(12) 


and 


<*u 


Vaj + o"y + a: — a v cr v — a v cr- — a v a- + 3 (t, 2 v + t + Ty r ) 


(13) 


Of the total strain energy under the stress-strain curve of figure 1 , the elastic strain energy com- 
ponent, represented by the triangle OBD , provides the static strength of the material. What makes the 
total strain energy important in the ultimate safety factor analysis is that local structural regions, having 
stress concentrations exceeding the elastic limit, will flow and deform to allow peak stresses to be shared 
by neighboring elastic regions. This benefit becomes obvious when comparing the limited elastic strain 
OD with the available plastic strain DF. Analysis must assure that local flow does not lead to global 
yielding or fracture. 


D. Triaxial Stress-Strain 

The inelastic stress-strain relationship may be developed from Hooke’s law, 

C| = - 4 ~ [CT, — |JL„(<T2 + CT 3 )1 

E a 

e 2 - ~r~ [02-M<*i+<*3)] (14) 

e 3 = -p- [03 +<T >)1 , 

E a 
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by replacing elastic constants with inelastic variables derived from uniaxial tests. E A is the inelastic 
modulus, and Poisson’s ratio, *x„, is stress dependent; this ratio is different about each principal axis, thus 
constituting an orthotropic material. The most practical approach is to let the Poisson’s ratio related to the 
largest stress represent the effective Poisson’s ratio for all axes. Other options and consequences are dis- 
cussed in section II-F. Assume cri is the largest principal stress, and *x„ = *i|. Squaring, adding equation 
(14), and collecting terms gives the combined stress-strain relationship over the entire inelastic range 


[£ i + £2 + £3] 


Ik. 

0+2*0 


= ai + ct 2 ‘ + ^3~ 


~ 4*i„ — 2*i„ ' 

1 + 2p.,i 


(cr,<T2 + 0’i0’3 + a 2 CT3) 


(15) 


The right side of equation (15) is the combined principal stress expression, and is noted to be 
similar to the distortion energy expression, equation (13), which is related to the uniaxial tension stress, 


<*A 



"I 1/2 


(<T|a 2 + 0-|CT3 +CT2CT3) 


06 ) 


Letting ct 2 = 03 = 0, then £ 2 = £3 = |i|£| and substituting into equation (15), provides the uniaxial 
relationships: 


(17) 

(18) 

Using equations (17) and (1), it follows that 



E a °a <t a 


e A 


p 

■57 “ Ea 


f e[±£i±ej 1 

L d+2*t„) J 


1/2 


or 


1 


E a 


K] (W 

[a) vb 


(19) 


Nondimensionalize equation (16) by dividing with af, and substituting into equation (19), gives 


1 


E a 



( 20 ) 
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where 


Cf) = 


1 + 


fflf + fcT- 4 > x „ - 2^,1 22 + Z1 + 

a,j L ai J L i+2 » jl « Jh a > ’ 


£2^3] 

O’? J 


\-b/2b 


( 21 ) 


But equation (21) was developed in its simplest form by assuming Poisson’s ratio independent of 
stress orientation, which is satisfied only at the elastic limit where b = 1 followed by <!> = 1 . Substitut- 
ing equation (20) into equation (14) yields the desired combined stress-strain relationships: 



The associated stress relationships, solved from equation (22), are 


£l 


£2 


£3 

<T| 


(a\' b ~E|(1 -p.|) + ^|(g 2 + S3) 

V'/ L 0 - 2 |ti)(l+M-t) . 


(*)"[- 


£ 2(1 ~ M-i) Mu(si T £ 3 > 

(1 -2fi,)(l +|i,,) 


] 


m 


£3(1 ~ M-l) + l JL l( e l + e 2) 

(1 — 2|JL,)( 1 + P-|) 


] 


(23) 


E. Shear Properties 

A shear stress-strain model is not only essential for analytical modeling, but for evaluating rosette 
strain gauges from structural test data. Published test data beyond the elastic limit is wanting. An 
approximation from uniaxial tension properties developed previously may be cautiously ventured. 
Recognizing that shear beyond the elastic limit is related to the same elastic-plastic slip phenomenon 
noted in figure 2, a two-parameter shear stress-strain expression is justified, 


t = c y d , 


(24) 
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and again the two parameters are determined from the extreme end conditions of the inelastic curve. The 
normal and shear stress-strain properties at the elastic limit end of figure 4 are established by 


Vel _ 
e el 


and 


Tel _ r _ 

Wl~ g - 


E 

2(1 + n) ’ 


(25) 


and related by the distortion energy equation (11) when all normal stresses are zero. 


T 




(26) 


Substituting equation (26) into the second part of equation (25) provides the elastic shear strain relation to 
the normal tension stress, 



d+1 x) 

E 


&EL - 


(27) 


The ultimate shear strain may be derived from the assumption that the ultimate shear strain energy 
is equal to the ultimate normal strain energy, 


r>lu 

J 7Tdy 


- J 


fit/ 


as 


(b+l) 


de 


Substituting equation (26) for the shear stress and integrating gives the ultimate shear strain required to 
determine the shear parameter in equation (24), 


7 u ~ e u 


yivi 

b + 2 


(28) 



Figure 4. Derived shear stress-strain model. 
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The two parameters are calculated from 


\0g((JulVEL) , v EL 

— — , and c = — — 

\og(yohEU V3 yd 


(29) 


F. Orthotropic Properties 

Polycrystalline materials experience orthotropic properties when multiaxial stresses exceed the 
elastic limit. It was noted earlier in figure 1 that stressing beyond point/!' and relaxing not only raised the 
elastic limit about one axis for the next cycle, but irrevocably used strain energy to cause permanent 
deformations about all axes. Properties about other axes remained the same as before uniaxial stressing, 
creating an orthotropic material. 

If a triaxially stressed article is tested beyond the elastic limit and relaxed, the article is essentially 
“sprung," and boundaries are permanently changed. The article will not behave exactly as it did before, 
when the external load cycle is repeating through the elastic operational range, and should be a considera- 
tion in prototype testing. If the prototype is disassembled and reassembled according to specifications, 
the prototype will behave elastically in the operating range as before testing. Beyond the elastic range, 
the fixed orthotropic properties are subject to change only when exceeded. 

In modeling inelastic substructures, boundaries are defined by the elastic global model. The 
substructure size selected must assure that remote inelastic stresses acting on the critical zone of interest 
do not spread to the elastic boundaries as the critical zone approaches the ultimate stress. 

Poisson's ratio is a property affected by plastic flow which is directional, and irreversible, and it 
increases with stress. In a triaxial stress field, a different Poisson's ratio is expected for each stress axis. 
The only restriction imposed on Poisson's ratio is that the volume growth. 


Ay 

Vo 


(l+e,)(l+e v )(l+e r )-l 


(30) 


at the elastic limit is established by the elastic Poisson’s ratio, and that the volume growth at fully plastic 
flow is zero. Modifying equation (22) to assume orthotropic Poisson’s ratio and substituting into equation 
(30) gives 


AV = 
Vo 



after dropping high order terms. 




(31) 


At the elastic limit, all Poisson’s ratios are equal to the elastic constant, and the volume growth 
becomes 
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(32) 


AV 

V„ 



(1-2M 


(g t . + g v + g.) 
(Tv 


At the fully plastic extreme, equation (31) is satisfied only if 


p, v = (x v = (UL r = 0,5 . 


(33) 


But these extreme limitations in triaxial stress are identical to those discussed for uniaxial stress, which 
implies that Poisson's ratios are independent of each other and dependent only on stress or strain along 
their respective axes. 

Since stress is the unknown variable in most structural problems, orthotropic Poisson's ratio 
related to multiaxial stress may sometimes complicate the solution. Generally, Poisson's ratio is not a 
very sensitive property, and a more practical approach is to assume the Poisson’s ratio related to the 
greater stress is common to all axes. If triaxial stresses are close, then a common Poisson's ratio will 
provide very good results. If stress magnitudes span the total plastic range, then Poisson’s ratios should 
be related to each respective strain where possible. Referring Poisson’s ratio to the average of two strains 
is another option. Sensitivity analyses of Poisson’s ratio as applied to critical models will identify the 
optimum option. 

Sections V and VI are examples of orthotropic properties applied to biaxial stress problems. Note 
that Poisson’s ratio in the denominator of equation (109) is related to dominant strain, e x , as suggested 
above. The maximum stress error that might develop from a maximum spread between e r and e- was 
estimated from the sensitivity equation (113). A 40-percent change in Poisson’s ratio (75-percent change 
in strain) caused less than 4-percent error in stress conversion. 


III. ANALYTICAL MODELING 


The basic analytical technique is similar to the computational piece-wise linear method often used 
to solve nonlinear problems. Unique to this analytical development was relating all nonlinear material 
models to engineering stresses and strains in exponential form, and synthesizing them with existing 
strength of material stress and strain models. This simplified the algebra and allowed continuous integra- 
tion for more descriptive results. The method uses existing principles of materials and strength familiar to 
most analysts, while circumventing plastic theory. It provides detailed insights to inelastic behavior. 
Models derived are related to the skirt structure and demonstrate the technique to be used for other 
models as required. 

The specimen of figure 5 represents a bar or plate having a length L and a rectangular cross sec- 
tion ol thickness t and width B normal to the page. External unit loads are the in-plane normal load, N , 
expressed in kips per inch, and the bending moment, M, expressed in kip-inches per inch. The bar speci- 
men model, whose width and thickness are relatively equal, is amenable to one-dimensional stresses and 
is a convenient introduction to the mechanics of inelastic behavior. 


13 




Figure 5. Specimen configuration and loading. 

A. Normal Tension and Bending Models 

Normal tension response models of rectangular bars is defined by uniaxial test parameters. The 
normal stress acting on the cross section of unit width is 



(34) 


and using equation (34), the corresponding strain is 



(35) 


Pure bending moment M is constant along the length L and deforms into a cylindrical surface of 
radius r shown in figure 6. Since cross-sectional planes have been shown as remaining plane in elastic 
and plastic bending, the centroid (t/2) of the rectangular cross section and the bending neutral axis are 
coincident from which the bending strain is linearly proportional to 



M \r^ 

Figure 6. Pure bending geometry. 
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Substituting tensile strain from equation (36) into equation (1) gives the bending stress distribution, 



along the cross section thickness as shown in figure 7. 


(37) 



The elemental area is dA = B dy and the moment resulting from stresses acting on the sectional 
area is dM = a y dA. Substituting the stress expression of equation (37) and integrating gives the total 
moment acting on the cross section. 



2aW2 )* + 2 
(b + 2)i* 


(38) 


The specimen width B is taken as unity, and hereafter omitted in most derivations. Solving for the curva- 
ture from equation (38) and then equating with (37) gives 


(b + 2) M 
2a(t/2) b+2 _ 



from which bending stress at the extreme fibers is reduced to 


(39) 




2(b + 2) M 


(40) 


15 





Note the elastic slope superimposed on the stress diagram, figure 7, is practically concurrent with the 
stress calculated from the nonlinear material model over the elastic range. This close fit is a measure of 
the accuracy expected from the analytical method. Also, note that when the strain hardening exponent, b , 
equals unity, the bending stress degenerates to the familiar elastic stress of a rectangular beam, as it 
should. 


Bending strains at the extreme fibers were similarly derived by substituting the strain relationship 
of equation (36) into equation (39), 


e m — 


2(b + 2) M 

at 2 




Mb 


(41) 


Bending transverse displacements may be derived from the moment and curvature relationship 
equation and from the slope expression, 0 = dyidx. The slope is related to the curvature by 


j_ = d% = < jy 

r dL dx 1 


Equating (39) with (42) gives 


d 2 y _ 

~d^~ 


] Mb 


and multiplying both sides by dx and integrating yields the slope expression of the midplane at x = 0, 

Mb 


0 = L r (6+2) m 

\2a(tl2f b + 2) 


(43) 


Integrating again gives the transverse deflection at x = 0, 


y = 


L?\ (b + 2) M 1 

2 \_2a(t/2) ib + 2) J 


Mb 


(44) 


B. Combined Stresses 

Shear-bending models have a wide variety of application and are essential for test analysis. Since 
no instrumentation is available to detect directly transverse shear V , a series of strain gauges measuring 
the variation of bending moment along the bar length L, figure 8, will provide the necessary data to 
determine the shear indirectly. 
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The moment, M, on the left interface is assumed to increase with distance, „v, to M + dM on the 
right interface. Induced stresses likewise increase from left to right. The stress-moment relationship is 
obtained by rearranging equation (37), 


_ a_ 

/ /> ’ 

and substituting into equation (38) to solve for the moment induced stress. 


( b + 2) / M 

2(t/2) b + 2 


Since ct v — ct m , the force acting on the left interface is 


o’ , dA 


(b + 2) / M 

2(tI2) b + 1 


dA , 


(45) 


(46) 


(47) 


and 


(a, + dcr x ) dA = - fc + 2 ^ > ’ (M + dM) dA , (48) 

2(t/2) b 2 

give the force on the right interface acting in the opposite direction. The net force is the difference of 
these two equations and is equal to the horizontal shear force in the m' — n' plane, 


T rv B dx 



(49) 
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Substituting the beam shear force V = dM/d.x. and the area dA = B dy into equation (49) and then 
integrating gives the flexural shear stress expression. 


Try 


V(fr+ 2) 

2 (b+ \)(t/2) b+2 


[(t/2) b+] -y b+i ] . 


(50) 


The horizontal shear is equal to the vertical shear along the cross section thickness on the right interface. 
Since the maximum shear stress of equation (50) occurs at the neutral axis (y = 0), where the bending 
stress is minimum, the maximum combined stress must be calculated from the distortion energy theory of 
equation (11). Shear stresses vary along the cross section, but unlike pure bending, the associated strains 
cause plane sections to warp, and should be considered when combining with other applied loads. 

Normal tension-bending stresses cannot be directly superimposed because they are not linearly 
related to strains, equation (34); nor is the bending neutral axis expected to coincide with the section 
centroid. However, an interacting model may be formulated by noting that normal strain e N is uniformly 
linear along the cross section, equation (35). Also, the bending strain is linearly proportional to the strain 
at the extreme fiber, s M , equation (36), because plane strains remain plane after pure bending. Since both 
axial strains are linear, they are algebraically added as shown in figure 9(A). 

+y +y 



(A) Strain Distribution (B) Stress Distribution 

Figure 9. Combined tension-bending distribution. 

The net strain distribution is linear and proportional to the intersection of both strains at a point C 
below the cross sectional centroid. The net strain is defined by 


e 



(51) 
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where e m is the maximum bending strain at the extreme fiber having the same polarity as the normal load 
N. The neutral axis, C, is negative in the y-direction. The half-thickness, t/2, is replaced by h for sim- 
plicity of notation, and the normal strain of equation (35) is calculated from the given applied normal 
load, geometry, and material properties. 

The combined stress distribution is shown on figure 9(B) and is determined by substituting equa- 
tion (51) into equation (34) to give 


a 


a 


£/V + 


(y-C) 

(h-C) 



sgn K^ 



(52) 


Expressions in absolute form allow raising strains to odd powers. SGN( ) is the signum function which 
reestablishes the sign of the expression. If its sign is positive, then the function equals + 1 and the strain 
is positive. If the function equals — 1 , the strain is negative. 

Solutions of the stress and strain distribution of equations (51) and (52) as a function of y rest with 
determining the unknown variables e M and C from the two static load equations; dN(y ) = ct(v) dy, and 
dM( y) = a(y) >’ dy. Substituting the stress expression of equation (52) into these static equations, the net 
normal load and total bending moment are 


N = af 


b ft 

dy + af 

p , (y-O 
e n + ,, e M 

- h 

(h-C) 

c 

(h-C) 


and 


II 

a 


b h 

y dy + a f 

.(y-C) 

E N + /t ^ e m 

-h 

(h-C) 

c 

(h-C) 


y d y 


(53) 


(54) 


respectively. Integrating within the limits, the normal load and moment equations become 


N = 


a(h — C) 
e M (b+ 1) 


Nil-t|»2] , 


(55) 


and 


where 


a(h-C) I 


~h+ {h '° 

~e N (/j+oi 

+ i!i1 


(b + 1 )bm 

IjiZ 

h (b + 2) 

_£m (/j-C)J 

I tu 1 

( b + 2) J 




(56) 


t|il = [AflS(e„ + £*,)]<* +,) , 


(57) 
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and 


4*2 


ABS 




(h± C) 
(A-C) 


- (/>+ 1 ) 


(58) 


Given the geometric property h and combined loads N and A/, the normal strain s N is calculated 
from equation (35). The bending strain at the extreme upper fiber e M and the bending neutral axis C are 
determined from equations (55) and (56), using the Newton-Raphson method presented in appendix A. 
Combined strain and stress distributions over the bar cross section may then be calculated by substituting 
the resulting bending strain and neutral axis values into equations (5 1 ) and (52), respectively. Conditions 
on figure 9 are that strains e N and e M have the same polarity, and that C is negative, but may be similarly 
modified for other conditions. 

Results from this analytical normal tension-bending model are shown in figure 10. The linear 
bending strain distribution exhibited in the strain diagram, figure 10(B), is superimposed on the uniform 
normal strain, as in an elastic beam diagram. However, the neutral and centroid axes are not coincident, 
and the location of the distance C between them is defined by equilibrium conditions in the nonlinear 
stress diagram, figure 10(C). There the areas of bending tension and compression representing loads are 
equal and opposite which cancel the axial load. Their moments about the neutral axis add to balance the 
externally applied moment. 


a 



Figure 10. Combined normal tension-bending diagram. 
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The unsymmetrical shape of the stress profile may be understood by projecting select points from 
strain and stress diagrams onto the uniaxial stress-strain curve, figure 10(A). These projected points are 
labeled PONM. It is clear that large tension strain from bending loads projected on the shallow slope 
between points N and M produces relatively small stresses, which accounts for the projected narrow ten- 
sion stress component into figure 10(C). On the other hand, small compression strains projected on the 
steep slope of NOP produce large stresses, which accounts for the very wide projected compression stress 
component. 

Stress and strain responses demonstrated in the diagram have interesting implications on struc- 
tural verification analysis. Back-to-back strain gauges are used to isolate pure normal strains from 
combined normal tension-bending strains. When testing in the elastic range, the pure normal strain is the 
average of two back-to-back strain gauge outputs. Beyond the elastic limit, this formula is not valid as 
revealed by the strain diagram. Though the bending stress is a third of the normal stress, the bending 
strain is over four times the normal strain. There is no direct method for isolating the normal strain from 
combined strain data except through an inelastic math model which is essential to support design, struc- 
tural test, and operational analysis. The combined normal-bending loads model provides a firm basis for 
predicting failure. It is a prerequisite for sensitivity analyses leading to robust designs and operational 
load placards to increase safety margins for unique missions. 


C. Plate Stresses and Strains 

All of the previous bar models were derived for rectangular cross sections whose specimen width 
dimension B along the Z-axis is not much greater than the thickness t. These models form a special class 
of one-dimensional stress specimen. Another type structure is the two-dimensional stress specimen, such 
as plates and shells, having length and width dimensions many times greater than the thickness, and in 
which the stress normal to the plate is considered negligible. A plain stress condition is assumed to exist 
in plate problems. 

Normal stress-strain relationships on plates may be derived by applying the plain stress con- 
ditions, o\. — 0, to the second and third of equations (23) to obtain the associated strain. 


e v = ~~~ (e v + £_-) . (59) 

1-p, 

Substituting equation (59) into the first and third of equations (23) gives the stress expressions in terms of 
biaxial strains. 


a, = a 




e, + p.e, "I 

l-M 2 J 


<j z = 


a(e, + |i,£ v ) 
[1-H 2 ]" 


fe v + |xe-] < ^- |> 


(60) 
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If the plate is subjected to uniaxial stress, in which a- = o\ = 0, then the first part of equation (60) 
reverts to e-/e v = — p, which defines the uniaxial Poisson’s ratio. The uniaxial strain is determined Irom 
the first of equations (22), 





which, again, reverts to the uniaxial expression of equation (1). It may be concluded that a plate in 
uniaxial normal stress behaves as a bar specimen in which the stress is delined as 


where h = t/2. 

Multiaxial property ( 1 — jx 2 ) was introduced into equation (60). Since the Poisson's ratio is related 
to the largest stress, it may be expressed as a two-parameter property for the compatibility reason given in 
section II, 


f v = (1-(X-) = Qa* • (62) 

Poisson’s ratio is calculated from equation (6) for end conditions \x,y at 07/ and p./,/ at o>/ , and substitut- 
ing into equation (62) gives 


V xU = [ 1 - M-vtyJ and ^ xhL - [ 1 - P-ix/J • 


(63) 


Both parameters of equation (62) are calculated as in section II, 


log 

t xU it xEI ._ 


log 

uS 

> 

b 

"2> 

l£j 



and 


Q = 


Tv/:/. 

OxF.L 


(64) 


Pure bending of plates results in a biaxial stress condition because the top surface is subjected to 
tension while the bottom is in compression simultaneously, and as in the bar specimen, produce lateral 
contraction and expansion from Poisson’s ratio effects, figure 1 1. However, plate widths are very large 
compared to the thickness, forcing lateral strains at the top, bottom, and middle surfaces to be equally 
restrained to e, = 0 for continuity [7]. 


Invoking the plate condition and bending constraint, a v = e_, = 0, into the third of equations 

( 22 ), 
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+<%, £ =0 



Figure I I. Poisson’s ratio effect on plates in bending. 


e, = 0 = — 


a v |(t- I 

M L^ _f T 


the plate lateral strain is satisfied when 


<T- 

^ = P-v 


(T 


(65) 


Substituting equation (65) into the first and second of equations (22) gives the other two plate normal 
strains 




l/b 


[ 1 - M-.v] , 


■ m 


(66) 


I lb 


[ — M-.v — M-.v 1 


Substituting equation (62) into the first of equations (66), the axial strain is rewritten as 


o = Q cr.r 


r H lb 

O , 


= Q 


W - 


(67) 


where the combined exponent is defined by 


1 +bR 


( 68 ) 
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Equating (36) with (67) gives the axial stress through the thickness 


a x 



P 



a 

r 


Q\ 


( 69 ) 


This stress acting over the plate cross section produces a moment, dM - cr v y dy. Integrating over the 
thickness, solving for the curvature, 




r l 

Gf 

$/b 

2_ = 

M 

0 + 2 

r 

2/, ( P + 2 j 

a 



i/p 


(70) 


and substituting equation (70) into equation (69) yields the plate bending stresses at the extreme fibers 


— 



(71) 


Substituting equation (70) into (36) gives the bending strain over the thickness. 



(72) 


and at the extreme fibers. 


— 


-h 

M 

’0 + 2" 

& 

$ib 


2h 2 

a 




(73) 


The slope expression at x = 0 is derived by equating (70) with (42) and integrating. 


0 = 


L 


M 


0 + 2 

2 h' p+2) 




Wp 


(74) 


Integrating again over the length, gives the deflection at x — 0, 

wu 

(75) 

The bending shear stress expression is analogous to the bar stress, using plate combined 
exponent, 0, 




M 

'0 + 2 ' 


0/6 


2/t ,p+J) 

a 
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( 76 ) 


v (ft+2)[ft (p+l) -y p+n ] 

20+1) h (p+2 ' 


Combined bending and normal loads interaction equations were derived similarly to the bar speci- 
men above. Substituting equation (61) into (1) gives the normal strain as a function of the known applied 
axial load, 


e n — 



(77) 


The combined normal and bending strains acting over the cross section are expressed by equation (51). 
The bending stress over the section is obtained from equations (36) and (69), 


— 


a 


( Mb) 


~\ 0 


Q 


ABS 


e M 


(y-c) 

(h-c) 


0 

SGN (y-c ) . 


The normal uniaxial stress is defined by equation (1), 


(78) 


c T/v = ae N h , 


which may be modified to conform with equation (78) by multiplying and dividing material parameters. 


CT/v — 


e N bl * Q 



0 


But 


e N b/& Q a< b -P> /b P 




so that the combined stress equation is simplified as 


0 - 7 ' 



ABS 


e A/+ e M 


(y-c) 


(h-c) 




(79) 


The net force acting on the cross section is dN(y) = fr-jiy) dy, and integrating over the limits — h 
to C, and C to +h gives 
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( 80 ) 


(h-Q [vp3 - 4*4] 
(0 + 1) e M 



M - 


p 1 

iJ/4 

1 

U h ~ C) 

(h + C) 

( 0 + 1 ) Em 

Q 


(0 + 2 ) 

Em ( h — C) 


+ »Ji3 



^+ ill 

(0 + 2 ) 



(81) 


where 



4/3 = 

[ABS (e„ + £ M )] ( 0 + '> , 

(82) 

and 

4/4 = 

/ (A+cnT ^ 0 

(83) 


V ( h-C)J 



Given the geometry, material, and applied axial normal and bending loads, the bending neutral 
axis C and bending strain e M are solved simultaneously from equations (80) and (81) through Newton- 
Raphson method presented in appendix A. The stress and strain profiles over the plate thickness are 
calculated from equations (51) and (79). 


IV. FEM VERSUS ANALYTICAL MODELS 


Two independent global and substructural models of the SRB aft skirt were developed using the 
ANSYS and NASTRAN commercial codes. Both substructural predictions of the critical weld region 
were in good agreement with gauge measurements on the structural test article throughout the elastic 
range. Beyond the elastic limit, the weld failed very short of predictions. Given that analytical methods 
solve few simple problems accurately and that FEM solves many complex problems approximately, the 
focus of the following analysis was to assess the nature and reasonableness of FEM approximations. 

The ANSYS finite element method was selected [6] because of its built-in graphic and contour 
mapping capability. The material behavior option used was the multilinear kinematic hardening (KH), 
and the element was a three-dimensional utilization brick (eight-node isoparametric) that was applied in 
the aft skirt structural model. Because a brick element is often used to fine tune predictions of triaxial 
stress zones having the highest potential for failure, it was the most critical and logical element to be 
explored. 

A common materials model was defined to be compatible with the ANSYS and analytical 
methods using equations (2) and (3), 


a = 87 e 


0.2 


(84) 
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based on assumed uniaxial properties <j f ± — 26.3 ksi, ct { ; — 45.7 ksi, - 0.04, and E — 10,500 ksi. 
The ANSYS code defines Poisson’s ratio similar to equation (6), 




0.5 — 0.2 


where E s is the secant modulus. The analytical method used the same formula which reduced to 


(a = 0.5- 95,000 ct 4 0 . 


(85) 


All other properties required for this method comparison are defined above and may be generated by the 
plate specimen program in appendix A. 

Two simple structures, as in figure 5, were modeled, both specimens having a length of 6 in and a 
thickness of t — 1 .44 in. The bar specimen had a width B of 1 in. The plate specimen had a width of 7 in 
to exercise the continuity condition in bending. Pure bending and tension loads were applied independ- 
ently and combined on both specimens providing sufficient number and variety of cases for comparison 
and cross checking. Loads were selected to produce stresses on the outermost fiber in excess of the elastic 
limit. 


A. Bar Specimen 

Results from the analytical method were based on equations (34), (35), (40), (41), (51), and (52) 
combined with the stress program in appendix A. Results from the FEM were based on data taken at the 
specimen midwidth and uppermost brick to represent typical stresses away from edge conditions. A 
comparison of the two sets of results, table 1 . shows that Poisson’s ratios compared well, but analytical 
bending strains ranged up to 20 percent higher than ANSYS. 


Table 1. Analytical versus ANSYS results of bar specimen. 


A 

B 

C 

D 

E 

F 

G 

H 

N 

M 

Ana 

lytical Me 
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Al 

MSYS Co 

de 

kip/in 

kip-in 
per in 

Stress 

ksi 

Strain 

in/in 


Stress 

ksi 

Strain 

in/in 


40 

0 

27.7 

.0033 

.34 

27.2 

.0032 

.34 

50 

0 

34.7 

.0101 

.43 

35.1 

.0100 

.44 

60 

0 

41.7 

.0252 

.47 

42.4 

.0228 

.48 

0 

15 

31.8 

.0065 

.41 

31.7 

.0062 

.39 

0 

18 

38.2 

.0163 

.45 

37.1 

.0139 

.44 

0 

20 

42.4 

.0276 

.47 

41.0 

.0225 

.45 

40 

2.5 

33.2 

.0080 

.42 

31.8 

.0066 

.41 

40 

7.5 

39.7 

.0197 

.46 

39.1 

.0182 

.47 

40 

10 

43.3 

.0304 

.47 

42.6 

.0275 

.48 
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Recalling that FEM data output relates to the average value of the outermost element, FEM bend- 
ing strain predictions are expected to fall short of surface prediction (fig. 12). This error increases with 
increasing moment and with decreasing number of bricks modeling the thickness. It is a characteristic 
error that applies to elastic bending stresses and strains as well as plastic strain computations. By 
extending the ANSYS linear strain plot to the surface, the modified strains came to within 5-percent 
agreement. This match on all three loading models, and notably on the intricate combined loading model, 
established a mutual credibility of the analytical and FEM approaches with a one-dimensional stress 
specimen. 


Characteristic 



B. Plate Specimen 

The ANSYS computational method was applied to the 7-in width specimen with no other change 
from the previous bar specimen computation. The analytical approach was similar to that of the bar speci- 
men, except that plate equations (77) through (8 1 ), combined with the plate program in appendix A, were 
imposed on the specimen. Results are presented in table 2. 


Table 2. Analytical versus ANSYS results of plate specimen. 
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40 

0 

27.7 
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.34 
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.32 

50 

0 

34.7 
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.43 

36.4 

.0089 

.35 

60 

0 

41.7 

.0252 

.47 

45.0 

.0204 

.29 

0 

15 

31.8 

.0059 

0 

32.2 

.0056 

.37 

0 

18 

38.2 

.0141 

0 

37.7 

.0105 

.41 

0 

20 

42.4 

.0231 

0 

41.9 

.0154 

.42 

40 

2.5 

33.2 

.0069 

.34 

32.5 

.0057 

.34 

40 

7.5 

39.7 

.0162 

.34 

40.7 

.0137 

.35 

40 

10 

43.3 

.0245 

.34 

44.3 

.0193 

.34 
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In comparing the plate specimen results, the FEM responses in columns F, G, and H due to axial 
normal loads only (M = 0) should be identical to those in table 1, since both specimens act in uniaxial 
tension, as discussed in section III-C. Contrarily, plate model strains in column G are shown lower than 
column D, and table 1 , implying a lateral constraint, e z = 0, where none was intended. If such a lateral 
constraint was imposed, strain ratios in column H would be naught, but in fact, FEM strain ratios are 
large. The ANSYS results of the uniaxially stressed plate as modeled are bewildered. 

The most naked departure between the analytical and FEM plate results was noted in the pure 
bending model for N = 0 and M = 15 through 20. It is generally established [7] that elastic plates in 
bending must assume zero lateral strain, e- = 0, at both surfaces and midplane away from edges in order 
to satisfy the continuity condition illustrated in figure 1 1 . The FEM lateral strains derived from column H 
are about a third of the axial strains, which clearly fails to comply with plate lateral constraint. This 
ANSYS fault is further manifested by the very small axial strain prediction in column G. Strains in 
column D represent the lower limits. 

The same surface characteristic error of figure 12 was observed in the FEM plate bending model 
which resulted in submarginal predictions. The standard remedy is to increase the brick elements across 
the thickness until the strains in surface elements converge. Since doubling the brick numbers more than 
doubles the machine time, subdividing only the surface elements within the aspect ratio limit should be a 
more prudent course. 

The net effects of plate continuity condition and FEM surface error may cause strain predictions 
to be 50 percent less than expected when compared with the analytical method. The ANSYS contour 
mapping feature does locate peak stress areas very well. Most patterns and trends may be valid, but 
intensities should be checked by other methods. 

Based on deficiencies demonstrated by the results of two different loading conditions, it seems 
that the ANSYS brick element as presently modeled with the plastic option may be inadequate to predict 
plate strains and associated stresses in critical two- and three-dimensional stress regions. 

The analytical method served to expose a conflict with the FEM plastic plate option, but the 
ANSYS elastic option with piece-wise nonlinear material input might have served as well. In fact, the 
ANSYS elastic piece- wise approach is recommended for primary inelastic computations (or backup) 
until the plastic option deficiency is resolved. Element number and sizes must be checked on all FEM for 
convergence of bending strains within reasonable tolerances. 


V. BUTT-WELD STRESSES 


Butt-welds are preferred joints for permanent assemblies because their mechanical behavior is 
uniform and undistinguishable from the connected base plate under normal operating loads. Beyond the 
elastic limit, dissimilar joint material properties and geometries promote stresses which influence the 
ultimate safety factor. The material discontinuity stress discovery was reported [ 1 ] in an exploration 
mode and is here synthesized into the ensuing analytical method. 
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Figure 13 is a square groove configuration of butt-welded plates. The plates are milled base 
material. The weld filler is usually the same base material with added substances to improve wetted 
surfaces and flow properties. The governing weld geometry is simply the thickness t limited by adjoining 
plate geometry requirements, and width w is usually decided by manufacturing considerations discussed 

before. 



Interfaces r 

Figure 13. Butt-weld specimen. 


A. Material Discontinuity 

Structural areas with greatest potential for yield and failure are often identified by severe, local 
discontinuities of loads, temperature, geometry, dissimilar base materials, and combinations. Another 
type of discontinuity is the insidious butt-weld of the same base materials; but interfacing materials 
exhibit sufficiently different mechanical properties beyond the elastic limit to cause local distortion. 
Property differences are primarily a result of the manufacturing process, such as HAZ gradient, casting 
semblance of filler, filler metal additives, work hardening, and heat treatment variations. In any case, the 
phenomenon and analysis at their interface are the same. 

When two such materials in tandem are subjected to a common axial stress, cr v , they exhibit a 
common strain and Poisson's ratio up to the elastic limit. Beyond the elastic limit, properties ot materials 
1 and 2 bifurcate to establish their independent stress-strain and stress-Poisson’s ratio relationships, 
which produces two different axial strains (E| and ct) <ind two Poisson s ratios (p,| and p. 2 ) as illustrated 
in figure 14. Because this difference of properties is the source of interface discontinuities, their relation- 
ship to stress must be redefined with their new origin at the bifurcation. 



Strain Poisson’s Ratio 

Figure 14. Properties of two materials under common stress. 
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The uniaxial stress-strain expression should be similar to the two-parameter form of equation ( 1 ), 


Gf.l.z) — ^f e .v — e EL.l\ B ' 


( 86 ) 


where subscript EL , 2 refers to the new origin, which is the elastic limit of weakest element, 2. Given 
ultimate and elastic limit stresses and strains of each element, and using equation ( 1 ) to calculate 
immediate stresses, o* j and a* >2 , at a suitable strain of weakest element. 


e*, 2 


e {;,2 

3 


(87) 


the exponent and coefficient of each element is computed as before. 


logUot; ~ CT/r L ,2 V(«* ~ «/■/ .2)] 
l0gl( e l/“ e /a.2V ( e k,2-£F.L ,2)] 


( 88 ) 
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[ £ a.2 — 


(89) 


The two-parameter stress Poisson’s ratio relationship is similarly developed. Using equations (5) 
through (9) and (1), Poisson’s ratios. 


* = *[C ■ <90) 

are calculated at ultimate stress for each element, p .u \ and p^ 2 - Intermediate stresses, oy* 1 and ay* 2 , are 
calculated for p,* = 0.36 from equation (90), 


°V* 



b/f 


The stress-Poisson’s ratio property beyond the elastic limit is 


( 91 ) 


(P-v P-*;/.) — (7 [a, a /v f , 


(92) 


and the two parameters are determined from 
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(93) 
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(94) 


The specimen of figure 13 is in a plane stress condition defined by o r = 0. Substituting this 
condition into the second of equations (22), the uniform lateral strain on either element subjected to axial 
stress, ct v , is 


e 



and after exceeding the elastic limit, the lateral strain assumes properties of equations (86) and (92), 


(e v — Bel.2) ~ ~C[ cr A — 0£x.2] / 


o\ 0 fL . : 


\/B 


which reduces to 
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(95) 


Lateral contraction on each element resulting from the lateral strain of equation (95) is calculated 

from 


n 

v = / (e y -e EL ' 2 ) dy 


(96) 


Applying equation (96) to the two elements of different materials subjected to a common axial 
stress, a lateral mismatch occurs at the free interface as shown, in figure 15(A) and expressed bv 


v, — v 2 = 
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Dividing by the last term, and neglecting the stress with the exponent having small differences between 
material properties of the same metallic base, the net lateral displacement of the two different materials 
elements at the free interface reduces to 


V(r = V | — v 2 = h 


G 2 


A 2< 


\/B2 


j [°a - Vel.i\ 


( F2XB2+ I VB2 


Cl A2 (l/fi2> 
G2 A l (,/fl,) 


(97) 


But elements I and 2 were bonded at the interface in an unloaded condition. In applying the axial 
load, the constrained interface caused local distortion as shown in figure 15(B). This distortion is the 
symptom of discontinuity stresses. A math model was devised in appendix B which consisted of an 
externally applied unit load q abruptly ending at a distance d, representing the interface of one material 
element. The abrupt load change simulated normal and shear discontinuity stresses at the boundary, 
which propagated through the weld (weakest) element. These discontinuity stress distributions were 
derived from classical mechanics, and results are stated in equation (B-9) in nondimensional form. 

The unit loading q of figure B-l acts normal (v-axis) to the element, and because of the 
orthotropic nature of polycrystalline materials, lateral material properties are independent of equations 
(86) and (92). Applying plane stress conditions to the second of equations (22), the lateral strain is 



(98) 


Substituting these parameters, and the first and second of equations (B-9) into (98) gives the 
displacement 




o 


[Dy - p-D.v] dy . 


(99) 


In reviewing the stress distributions along the half-thickness in figure B-3, the D x component is 
noted to describe equal positive and negative areas, which implies that the intergration over the thickness 
would negate its contribution to lateral displacement in equation (99). Deleting D x and integrating gives 
the lateral displacement at the boundary as a function of the dummy loading. 
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where 
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Figure 15. Material discontinuity displacements. 


To match the net displacement of figure 15(b), the external dummy unit load in equation (100) 
must compress element 1, and stretch element 2 with equal and opposite intensity. 


v — V] + Vi — 1 H — — — 

q " a <•/«,) Ks a 2 


Substituting equations (97) and (102) into the bonded interface equilibrium condition, v<r + v q 
0, the external dummy-to-axial unit loads ratio is 


34 



( 103 ) 
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Finally, multiplying equations (103) and (B-9) eliminates the dummy load q and produces the 
desired discontinuity stress distributions at the interface as a function of applied axial stress beyond the 


elastic limit. 
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where a, > a> 7 2 and a> 7 2 refer to the weakest element, 2. Most relevant to a designer is the maximum 
stress factor occurring at the applied ultimate stress of the weakest material. Equations (107) are the 
ratios of stress components with the ultimate stress: 
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Figure 16 presents the discontinuity stress factors distribution along the interface from mid- 
thickness to surface, as defined by equation (107), using properties from figure 14. The weld geometry 
assumed was that of the SRB skirt, t/w = 5.75. Note that the shear stress component is largest and occurs 
at y/h = 0.96. It is not possible to experimentally measure this shear because it is transverse, and because 
it occurs just below the surface. The shear distribution is very sharp, and it required 42 points to focus the 
peak and would have eluded most FEM modelers using six element thickness selection. The peak value 
was compromised further because an element averages the stresses over each brick thickness. Never- 
theless, shear presence was hypothesized by Mr. Paul Munafo from a metallurgical viewpoint during the 
skirt failure review, but it could not be supported by mechanics at that time. Because material models are 
so arduous, a program to calculate them and stresses of equation (107) is presented in appendix B. 

The failure expression of equation (11) may be used to compute the combined discontinuity stress 
factor at the interface. Including a factor of one onto the .v-stress component to account for the externally 
applied axial stress, the combined stress factor is 


>Uc 


-VI 


1 + S(Jx) 2 + Si/y — ( 1 + S(Jy) S(Jy + 3 S {Jyy . 
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The material discontinuity stress factor for the above butt-weld example was 1.11, which 
occurred 1 /32-in below the weld surface. 



Location Above Mid-Plane, y/h 

Figure 16. Discontinuity stress factors along weld interface. 


B. Geometric Effects 

Benefits of improved specific strength materials in high performance shell structures are often 
tempered by joint weaknesses and further discounted by joint geometric limitations. Butt-weld simple 
geometric parameters were examined for inelastic mechanics effects using models developed in sections 
II and III. 

Applied Inplane Stresses: Discontinuity stresses were calculated for plate specimen of butt-welds 
having a thickness-to-width ratio of t/w = 5.75, figure 16. Other ratios were varied over a practical range 
to ascertain the existence of a weld geometry that would decrease the peak stress factors at the interface. 
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The governing parameter was a stress factor Kq of equation (104). Using the same material properties 
and location of maximum combined stress, y/h = 0.96 observed in figure 16, a combined discontinuity 
stress factor of 1 . 1 1 was noted to occur over a range of t/w = 1 to 7. Based on these results, there is no 
requirement from the mechanics viewpoint to control it. The overwhelming geometric consideration 
seems to be the heat sink effects and associated thermal distortions induced into the structure. 

Test specimens of weld joints are usually shaped with nearly square test cross sections which 
impose triaxial stress conditions and with unique material discontinuity stress factors at the weld 
interlace. However, most butt-weld joint applications are on plates having cross sections of very large 
aspect ratios which result in a two-dimensional stress field causing the discontinuity stresses defined 
above. From the designer's viewpoint, it is necessary to know that current test properties from square 
sections are not more benign than those experienced in plates and shells. Projects should incorporate a 
plan to test a statistical sample size of both configurations to provide the designer with an appropriate 
data base. 

Bending Stresses: In-plane loading is the preferred strength application of butt-welds, but these 
same loads can induce bending when welds are designed into regions with abrupt geometric changes. 
This principle is especially appreciated in girth welds of pressure vessels where abrupt changes in 
meridional curvatures between end-closures and cylinders cause abrupt displacement differential at the 
intersection. Displacements are matched through shears and moments which induce local discontinuity 
stress waves that damp out over a length related to ( R t) n . In particular, local bending causes materials 
with the lowest elastic limit to hinge and assume a disproportionate share of distortion. The combination 
of the peak stress waves falling in the girth weld having the lowest elastic limit causes the weld material 
to yield first and progressively distort most. 

Because bending displacement to failure is a function of weld width noted by equations (74) and 
(75), the wider the weld or the HAZ, the greater the safety margin. But controlling this geometric 
parameter is not a reliable method for designing high performance structures. Consequently, good design 
practices suggest that all welds be designed away from discontinuities wherever possible. 

If local bending cannot be designed out of butt-welded plates, then the weld width determined by 
manufacturing process discretions must be analyzed to assure sustenance of the superimposed local 
moment displacements to within the specified margin; if it cannot, the plate thickness should be 
increased. Nevertheless, the final decision on weld width should invariable continue to rest with metal- 
lurgists and manufacturing considerations. 


VI. STRAIN GAUGE DATA ANALYSIS 


Strain gauge data analysis closes the loop on the structural modeling process. Structural failure 
predictions are based on four critical design uncertainties; maximum expected operating load analysis, 
manufacturing and processing quality and tolerances, strength of material confidence, and stress 
modeling assumptions and techniques. Only the last three combined uncertainties are verifiable from 
development structural test through strain response measurements on a test article under well-defined 
applied loads. Preferred measuring instruments are electrical strain gauges. 
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The primary purpose for applying the multitude of gauges that often move project managers is to 
assess the structural behavior of critical regions, to correlate measurements with the stress model predic- 
tions, and to ultimately fine-tune the prediction model. The calibrated model is then used to extract and 
verify loads, again through strain gauge data under actual operational environments. At this point, the 
uncertainties are not only resolved, but the verified model is used to determine safety margins with a 
variety of projected operational loads. 

Another important application of strain gauge data correlation with prediction models is the 
verification of prototype safety factors. Because prototype test structures must remain operational after 
test, a maximum test load of 20 percent above the elastic limit is allowed. In this 20-percent nonlinear 
response sample, the strain gauge data is tracked with the nonlinear prediction model. If the gauge data 
and model predictions correlate very well, then the prototype article is expected to support the predicted 
ultimate safety factor. 

In all the above cases, strain gauge data evaluation is a direct feedback to modeling predictions. 
To begin, strain gauge types, ranges, locations, and orientations are determined from a map of critical 
stress regions based on prediction models. The same map is required for evaluating the resulting test 
strain gauge data. Recorded strains are converted to the conventional domain of stresses used specifically 
in prediction and materials data bases. The conversion of multiaxial strain to inelastic stress variables is 
made difficult because the Poisson’s ratio increases with stress, which in turn is nonlinearly related to 
stress. A conversion method follows which is based on the preceding analytical approach. 

Stresses normal to a plate, or shell, are considered negligible so that a plain stress condition may 
be assumed to exist, o\. = 0. The related principal stresses are similar to those derived in equation (60), 
and rewritten with orthotropic Poisson’s ratios. 
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Poisson’s ratios are related to strain about each axis by equation (5), 
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Substituting equation (111) into equations (109) and (110) produces the required principal stress 
equations expressed in measured strains. 
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Solving for the dominant stress, a v , from the first part of equation (112) and substituting into the second 
equation provides the two principal normal stresses. 

Note that Poisson’s ratio in the denominator of equations (109) is related to dominant strain, e v , 
as suggested in section II-F, though it is shared by strains e, and e-. The maximum error caused by this 
simplification may be determined from the nonlinear expression of Poisson’s ratio in the first of 
equations (112). Differentiating stress a v with respect to e, in the denominator and dividing by the stress 
and strain gives the percent error of stress with percent change in strain. 


da. r _ 200bfge 2/ de x (113) 

CT » l - g E y e x 

where the maximum ds x Je x = (e v — e.)/e v . 

Equation ( 1 12) is suitable only when axes of principal stresses on the test article are known and if 
shear stresses are not existent. A set of perpendicular strain gauges, tee rosettes, are used to measure 
these principal strains. Symmetrical pressure bottles and shells of revolution are examples of zero shear 
presence; but regions on plate or shell known to have inplane shear stresses must be instrumented with 
rectangular rosettes, figure 17, to obtain the third equation for solving three inplane strains. 



The rectangular rosette has two perpendicular strain gauges with a third bisecting them at 45 
degrees. They are mounted at some known angle with respect to a reference .v-axis marked on the test 
article as shown. The relationship between strain gauge readings and the normal and shear strains is 
defined by the transformation matrix, 
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where the Roman numeral subscripts refer to the measured strains. Equation (1 14) has the form 
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X = [A] Y , 


from which the reference axis strains are resolved through its inverse matrix form. 


Y = [A]-' X . 


( 115 ) 


Restating results of the stress tensor in section 11 with strains and substituting the strains resolved 
from equation (115) provide the principal strains at the rosette location 

S x + E z J . / T~ 2 T 7 2 

fi|,2 = - 2 -± ~2 Vfiv +e -~ - 2e ' £ -- + 7v_- • (116) 

The angle producing principal strains is 


tan 2 0 /; 


7xz_ 

e v — e. ’ 


(117) 


and the maximum shear strain is 


"Ymax = V e.v 2 + e . 2 - 2e t e- + y„ 2 . (118) 

To convert principal strains to principal stresses, substitute results from equation (116) into equation 
(112) as before. 


VII. MODELING VERIFICATION 


A. Materials 

There is adequate documentation on normal stress-strain data through the inelastic regime, but 
none on Poisson’s ratio and shear strain at ultimate stress to verify models in section II. Orthotropic 
Poisson’s ratios may be backed out of biaxial strain responses from known stresses developed on pres- 
surized elliptical domes. This same test could verify strain gauge conversion equations (112). 

Another task would be to verify the materials discontinuity stress phenomenon. A simple 
approach would be to test welded plate specimens having a width seven times the thickness to failure and 
compare the strength results with current weld specimen having square test sections. A more 
comprehensive task would incorporate a photoelastic coat to assess plastic strain variations associated 
with HAZ, fillet, and base properties concurrently in the weld region. 
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B. Analytical Models 


Though analytical model developments were based on logical extensions of sound physical 
principles, some simplifications were necessary to reduce certain models to practical forms. Except for 
the materials discontinuity stress development, a necessary mathematical check for all inelastic models 
was to substitute the linear parameters a = E and b — 1 to assure these expressions collapsed to the 
familiar textbook strength of materials formulas. It should be interesting to experimentally verify the 
plate bending model. 

Since SRB skirt ultimate safety factor prediction by the NASTRAN code was similar to the 
ANSYS code, using both plastic options and brick elements, it would seem appropriate to also suspect 
NASTRAN plate results until verified by a similar ANSYS task [6]. 


VIII. SUMMARY 


The safety of all structures depends upon the inelastic strength to ultimate stress behavior, yet this 
phenomenon is not well known. Perhaps there has been little incentive to dwell into intellectual nonlinear 
analyses that are bumped against arbitrary safety factors when a structural test to failure can more posi- 
tively confirm the ultimate safety margin. However, it does not provide the physics of failure. When a 
critical high performance structure failed short of predictions and a redesign again failed short with no 
further clues as to the mechanics of failure, it became necessary to understand the structural behavior of 
inelastic stress as applied to design through test analysis. 

The behavior of structures beyond the elastic limit is typical of the elastic-plastic mechanics of 
appendix C, which is neither very versatile nor comprehensible to the prevalent designer. Instead, the 
approach employed was an analytical method similar to the piece-wise linear technique used with linear 
computational codes and stepped nonlinear material inputs. The analytical method difference is that non- 
linear material properties were expressed in engineering stresses and strains and modeled in exponential 
form, and were then integrated with strength of materials linear stress and strain expressions. The inte- 
gration of these exponential expressions is continuous which allows for fine focused solutions. 

Several analytical structural models were developed that are basic to plate and shell welded con- 
struction, and predictions were compared with a commercial FEM. A common strain error characteristic 
to all FEM bending models was noted. Computed maximum bending strains of FEM codes are not at the 
surface, but they are predictions at the centroid of the surface element, figure 17, resulting in a lower 
strain prediction. 

A simple plate structure was modeled with the same ANSYS code using the same SRB skirt 
plastic material and brick options. The uniaxial specimen produced mixed results. Axial strain predic- 
tions were low enough to imply an inappropriate lateral constraint, but the strain ratios did not comply 
with nor without that constraint. Bending specimen results did not conform to plate continuity conditions 
in which the lateral strain along the width must equal zero. In comparing the net results of the ANSYS 
plate compliance and surface error predictions with the analytical method, the ANSYS method rendered 
50 percent less strain. 
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Based on the few SRB related examples, the ANSYS plastic brick option was clearly not an 
appropriate model for the aft skirt stress analysis beyond the elastic limit. Since the NASTRAN code 
predictions, using similar options, fell just as short as the ANSYS, the NASTRAN method should be 
suspect until scrutinized similarly to the ANSYS task [6], Inelastic analysis may still be accomplished 
with either code by incorporating nonlinear material properties into the linear elastic option in the 
familiar step-wise fashion. This step-wise method is recommended as a back-up approach for future 
design analysis. 

Discovery of weld discontinuity stresses was reported earlier [1] and was incorporated into the 
analytical method developed in this study for completeness. Results indicate shear stress is a significant 
factor in the failure criteria, which has a very narrow distribution, and peaks 2 percent below the plate 
surface at the material discontinuity interface. It is not possible to experimentally detect it, and it is dif- 
ficult to model it through FEM. Combined discontinuity stresses are not dependent on geometry, but are 
strongly dependent on the relative properties of the interfacing elements; first on the ratio of their 
elastic limit stresses, and then on the strain hardening exponent ratios. The implication is that material 
discontinuity stress may occur in a very narrow band of the weakest material. 

The exact stress and strain coupling under combined tension and bending loads was resolved 
through the analytical method with some rather surprising behavior. Figure 10. When measured strains 
from back-to-back gauge exceed the elastic limit, the bending strain may be several times larger than the 
axial strain, but the stress is a small fraction of the axial stress. This model provides the basis for assess- 
ing operational loads sensitivities. Because it is an often encountered condition and lengthy to solve, it 
was programmed and included in appendix A. 

Another surprise was the tee rosettes strain gauge conversion to stress expressions. While one of 
the stress formulas was similar to the elastic conversion, the second, equation (1 10), was distinctly dif- 
ferent. In substituting elastic properties into equations (109) and (1 10), both equations collapse to the 
familiar elastic conversion formulas. 

Analytical models presented are related to the skirt failure, other models may be similarly derived 
as required. What began as a troubleshooting task evolved into an analytical technique waiting to happen. 
To understand the inelastic behavior of specific structures, existing nonlinear material models were 
modified to combine with existing strength models. As more models were derived, a common technique 
emerged that was extendable to combined and multiaxial loading. When the technique was not only 
applicable, but required, to convert strains from test data to stresses, the analytical loop was complete. 
There are limits to its application, as there are to plasticity theory, but the technique is another option 
friendly to a larger spectrum of analysts. 


42 



REFERENCES 


1. Verderaime, V.: "Weld Stresses Beyond Elastic Limit, Material Discontinuity." NASA TP 2935, 
August 1989. 

2. Bianca, C.J.: Presentation on "Use of Safety Factor of 2.0 for No-Test Option." January 24, 1990. 

3. MIL-Handbook C5, “Metallic Materials and Elements for Aerospace Vehicle Structures," 
September 1976. 

4. Vaughn, T.P.: MSFC Metallic Materials Division Project 86-020-EH23, report to be published. 

5. Freudenthal, A M.: "The Inelastic Behavior of Engineering Materials and Structures." John Wiley 
and Sons, 1950. 

6. Berry, D.M.: “Aft Skirt Analysis." Vol. I, Contract NAS8-37282, April 1989. 

7. Timoshenko and Gere: “Theory of Plates and Shells.” McGraw-Hill, second edition. 

8. Fliigge, W.: “Handbook of Engineering Mechanics.” McGraw-Hill, 1962. 


43 




APPENDIX A 


NEWTON-RAPHSON METHOD 


1. Combined Stresses, Bar Specimen 

Two unknown variables. C and e u . are implicit in the two equations (55) and (56). These vari- 
ables may be determined through the Newton-Raphson method by rearranging both equations in the 
following functional form. 
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and estimating starting values for these unknowns. Then. 
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PROGRAM. COMBINED PLASTIC AXIAL-BENDING LOADS (BAR SPECIMEN) 
moment must be greater than zero 


' INPUTS 

M=7.5 

P=40 

H=.72 

A=98 

B=.22 

BB=B+1 

PH= P/(2*H* A) 
EP=(ABS(PH)) A (1/B) 


'external applied moment, kip-in/in 
' " axial load, kips/in 

' plate half-thickness 
'strength coeff. of mat'l 

■hardening exponent of mat! 

'strain caused by axial load 


'estimate unknown variables 

INPUT "estimate EM=";EM 'strain caused by moment 

INPUT "estimate C='*;C Trending neutral axis from y=0, inch. 0<= -C<-H 


' NEWTON-RAPHSON METHOD 
J=0 

START: Ql=EM*.0011: Q2=C*.0011 
CALL FUNCTION(M,P,H,A,B,BB,EP,EM,C,Fl ,F2) 
D1=EM: D2=C: Y1=F1: Y2=F2: EM=EM+Q1 
CALL FUNCTION(M J>.H, A,B ,B B JEP.EM.C J 1 ,F2) 
G1=(F1-Y1)/Q1: G2=(F2-Y2)/Q1: C=C+Q2: EM=D1 
CALL FUNCTION(M,P t H,A,B,BB»EP,EM,C,Fl,F2) 
L1=(F1-Y1)/Q2: L2=(F2-Y2)/Q2 


'cramer's rule 

K»G1*L2-G2*L1 

K1=Y1*L2-Y2*L1 

K2=G 1* Y2-G2*Y 1 

EM=D1-K1/K 

C-D2-K2/K 

J-J+l 
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PRINT "EM«",EM," C=\C 
IF ABS((D1-EM)/D1)>.000001 GOTO START 


'AXIAL- BEND PROFILE 
OPEN"CLIP:" FOR OUTPUT AS #3 

YS-H: YF=-H: NY=25 
N-NY-1 
DY=(YF-YS)/N 
ET=0: ST=0 
Y*YS 

FOR 1=1 TON 

ET=EP+(Y-C)*EM/(H-C) ’HT= combined strain 

ST=A*((ABS(ET)) A B)*SGN(ET) ’ST* combined stress 

SP=P/(2*H) SP = AXIAL STRESS 


WRITE #3,Y,ET,ST,EP,SP 
PRINT Y, ET.ST3P.SP 
Y*YS+(I+1)*DY 
NEXT I 

CLOSE #3 
STOP 

SUB FUNCTION (M,P,H.A,B,BB,EP,EM,C,F1,F2) STATIC 

R1=A*(H-C)/EM/BB 

R2-EP+EM 

R3=EP-EM*(H+C)/(H-C) 

R4=((ABS(R2)) A BB) 

R5=((ABS(R3)) A BB) 

R6*H+((H-C)/(B+2))*(EP/EM-(H+C)/(H-C)) 

R7=H-((H-C)/(B+2))*(EP/EM+1) 

FI »R1*(R4-R5)-P 
F2 »R1*(R4*R7+R5*R6)-M 

PRINT EM, C 

END SUB 


2. Combined Stresses, Plate Specimen 


The two unknown variables. C and e M , implicit in equations (80) and (81 ) may be solved by the 
Newton-Raphson method. Rearranging both equations in the functional form as before. 



proceed with equations (A-3) through (A-9) above. A program used in section IV in this document 
follows. 
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PROGRAM, COMBINED PLASTIC AXIAL-BENDING LOADS ( PLATE SPECIMEN) 


’ INPUTS 

M=10 

P=40 

H=.72 
SU=45.7 
SE=26.3 
EU=.04 
E= 10500 

B=LOG(SU/SE)/LOG(EU*E/SE) 

A=SE/((SE/E) A B) 

BE=1/B 

SP= P/(2* H) 


external applied moment, kip-in/in 
' " axial load, kips/in 

' plate half-thickness 
'max stress, ksi 
'elastic limit stress, ksi 
'max strain 

Young's modulus, ksi 

'hardening exponent of mat'l 
'strength coeff. of mat’l 

'applied axial stress, ksi 


"Poisson's ratio 
UUA»((SE/A) A BE)*E/SE 
UA=.5*UUA+.3*(1-UUA) 
TJ=.5-(UA-.3)*S/E*((A/S) A BE) 

'plate coefficient 

UXEL=.5-(UA-.3)*SE/E*((A/SE) A BE) 
UXU=.5-(UA-.3)*SU/E*((A/SU) A BE) 
G AEL= 1 -UXEL A 2 
GAU*1-UXU A 2 

R=LOG(GAU/GAEL)/LOG(SU/SE) 

Q=GAEL/(SE A R) 

BP=B/(1+B*R) 

EP=(SP/A) A BE 


'poisson's ratio at any stress s 

'poisson ratio at elastic limit 
'poisson ratio at ultimate 

'plate coef at elastic limit 
'plate coef at ultimate 

'plate combined expo. 

'strain caused by axial load 


'estimate unknown variables 
INPUT "estimate EM=";EM 
INPUT "estimate C=";C 


'strain caused by moment 
'bending neutral axis from y=0, inch. 0<= -C<-H 


' NEWTON-RAPHSON METHOD 
J=0 

START: Q1=EM*.0011: Q2=C* .0011 
CALL FUNCTI0N(MJP,H,A,B3E,EP£M,C,BP.Q JF1 ,F2) 
D1=EM: D2=C: Y1=F1: Y2=F2: EM=EM+Q1 
CALL FUNCTION(M.P,H,A,B,BE,EP,EM,C,BP,Q,Fl,F2) 
G1=(F1-Y1)/Q1: G2=(F2-Y2)/Q1: C=C+Q2: EM=D1 
CALL FUNCTION(M,P,H,A,B,BE,EP,EM,C,BP,Q,Fl,F2) 
L1=(F1-Y1)/Q2: L2=(F2-Y2)/Q2 


'cramer's rule 

K=G1*L2-G2*L1 

K1=*Y 1 *L2- Y2*L 1 

K2»G1*Y2-G2*Y1 

EM=D1-K1/K 

C-D2-K2/K 

J*J+1 


PRINT "EM=",EM," C=",C 
IF ABS((D1-EM)/D1)> .000001 GOTO START 


48 


•AXIAL- BEND PROFILE 
OPEN" CLIP:" FOR OUTPUT AS #3 



YS=H: YF=-H: NY=25 

N=NY-1 

DY=(YF-YS)/N 


ET=0: ST=0 
Y=YS 

FOR 1=1 TO N 

ET= EP+(Y-C)*EM/(H-C) 'total strain vs thickness 

ST=(((A A BE)/Q) A BP)*((ABS(ET)) A BP)*SGN(ET) 'total stress vs thickness 


WRITE #3,Y,ET,ST,EP,SP 
PRINT Y, ET,ST,EP,SP 
Y=YS+(I+1)*DY 
NEXT I 

CLOSE #3 

STOP 

SUB FUNCTION (M,P,H,A,B,BE,EP > EM,C,BP,Q .F1.F2) STATIC 
R 1 =(((A A BE)/Q) A BP)*(H-C)/EM/(BP+ 1 ) 

R2=EP+EM 

R3=EP-EM*(H+C)/(H-C) 

R4=((ABS(R2)) A (BP+ 1 )) 

R5=((ABS(R3)) A (BP+1)) 

R6=H+((H-C)/(BP+2))*(EP/EM-(H+C)/(H-C)) 

R7=H-((H-C)/(BP+2))*(EP/EM+1) 

FI =R1*(R4-R5)-P 
F2 =R1*(R4*R7+R5*R6)-M 

PRINT EM, C 

END SUB 
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APPENDIX B 


DISCONTINUITY STRESS MODEL 


A model is required for determining local stress components on a fillet weld specimen caused by 
an abrupt change in lateral displacement at the bonded interface. Induced inplane shear stresses. t vv . 
should be zero at the top and bottom surfaces and at the midthickness, and the lateral normal stress. ay. 
should be zero at both surfaces. These boundary conditions may be satisfied by inventing a dummy 
interfacing element having the same thickness as the fillet and subjecting it to a uniform external load. q. 
over a suitable length, d. The load, q, induces a uniform lateral strain. e v . across the thickness except at 
the interface, where the stresses and displacements are matched with the fillet element. 

Since the dummy element is used only to simulate boundary conditions at the weld bond, stress 
distributions for.v < d are not relevant to the problem. Therefore, the lateral loading, q. and axial load- 
ing. o\. will be related to the common interface displacement to eliminate the lateral load from the dis- 
continuity stress analysis in section V. It follows that the origin of the v-axis will be shifted in the text 
from the position shown in figure B-l to a point on the interface. 


y 



Figure B-l. Weld interface model. 

Figure B-l describes a simply connected contour with inplane forces which makes it independent 
of material constants. The geometry and the discontinuous dummy loading q are symmetrical about the .v- 
and y-axes. This is a classical problem [8] applicable to the generalized plane stress method. It consists of 
integrating the compatibility equation. 


8 4 <t> ! 2 d 4 ^ ( 3 4 <fr = 0 

dx 4 dx 2 dy 2 dy 4 


(B-l) 


where Airy's stress function. <t>, is defined by three components derived from equilibrium: 




a 2 d> 

dy 2 


a v 


a 2 d> _ a 2 4> 

dx 2 ’ <T " dxdy 


(B-2) 
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Equation (B-l) is of the form <t>(.v. v) = XY whose solution is well known to be 


0(.v,v) - (C| cos X.v + C 2 sin X.v) (C 3 cosh Xy + C 4 sinh Xy + C s v cosh Xy + C 6 y sinh Xy) (B-3) 


where Ci,...C 6 are arbitrary constants and X = rmr/L for m = 1,2,... Substituting equation (B-3) into 
(B-2), satisfying stress boundary conditions 

t vv = -- - - -- = 0 at x — ± h and v = ± L , 
dxdy 

and allowing that sin A_\ to vanish for a symmetrical loading of 

e = /t +is cos . (B-4) 

l L 

at y = ±h, determined the constants to be 


C i = I and C 2 = C 4 = C 5 = 0 . 


Stress components of equation (2) become: 

-2 i 

= — — = 2 Q S ^v ' , COS . V , [(X/i cosh X/t — sinh X/t) cosh Xv — (sinh X/?) Xy sinh Xy] , 
dy 2 2 X/t + sinh 2 X/t 

oy = -- y - = 2Q S cos [(sinh X/i) Xv sinh Xv — (X/i cosh X/t + sinh X/t) cosh Xy] , 

dx 2X/t + sinh 2X/t 

T rv = ~ f ^ ' = 2 (? 2 S * n rr . [(X/? cosh X/») sinh Xv — (sinh X/i) Xv cosh Xv]. (B-5) 

dxdy 2\h + sinh 2X/t 

Dummy uniform tension loading terms, A and B , of equation (B-4) are found in the usual way: 


f A dx = f q dx , 

-z. -a 

where A - qd/L , and 

B = y- f q cos X.v dx + f (0) = sin Xt/ . (B-6) 
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Evaluating A for integers 0 < m < I in each stress component as \ approaches zero for v = 0, A becomes 
zero for cr v and t vv , and AH for cr v . Let nn t = a and 


\ci — ad/L — a T , \x = o uc/L — aX , \h = ah/L = a H , Ky = a y/L = a Y . (B-7) 

Substituting the evaluated A, equation (B-6), and identities (B-7) into (B-5) gives the discon- 
tinuity stress factors in nondimensional form, 


— = 4^ — sin aT cos aX [(a.H cosh a H— sinh a H) cosh aK — (sinh aH) a Y sinh aK)] 

(2aA/ + sinh 2a H) a 

— = T+4 2 — sin aT cos ccX [(sinh a H) a Y sinh aY — (a.H cosh a H+ sinh aH) cosh aK)] 

q (2a// + sinh 2aH) a 


= 4 2 sin aT sin aX [(a// cosh aH) sinh aY - (sinh aH) aY cosh af)] . 

Q (2a// + sinh 2a H) a 

(B-8) 


Discontinuity stress factors distribution over the length and half thickness are calculated from 
equations (B-8) and plotted in figure B-2 for m = 4 which was limited by program overflow. Discon- 
tinuity stress factors distribution at the bond interface are determined from equation (B-8) by letting aX 
= aT. Recognizing that aH is large so that sinh aH = cosh aH, and 2a H is small, discontinuity stresses 
at the interface are 


m F/ \ 

D, = ^ = 4 S sin aT cos aT |e~ a// ] U/--1J cosh aY-Y sinh aE 


m p . \ 

D y — — = T+4 2 sin a T cos a T [c " W J Y sinh aK— (// + — jcosf 

q i V a / 


D xy — = 4^ sin a T sin a T [ e U "\[H sinh aK— Y cosh aK] 


(B-9) 


Stress component factor distributions of equation (B-9) are plotted in figure B-3. 
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PROGRAM MATERIAL DISCONTINUITY, FIGURE (16) 

STRESS FACTOR @ ULT STRESS DISTRIBUTION ALONG Y-AXIS INTERFACE 

Hl=h/1, H2=d/I, Yl=y/1, X2=x/1, SXY=shear, SX=axial. SY=norm. i/w = HI/ (1-d) 

OPEN "CLIP:" FOR OUTPUT AS #2 


' geo input 
H2= .2 
HI =4.6 
Ml= 5 


' material INPUTS 

SU1=58 

SU2=48.6 

SE1=32 

SE2=26.5 

EU1=.04 

EU2=.04 

E= 10500 

EE1=SE1/E 

EE2=SE2/E 


'max stress PLATE, ksi 
DO WELD 

'elastic limit stress PLATE, ksi 
DO WELD 
'max PLATE strain 
' DO WELD 


Young's modulus, ksi 
'elastic limit 


' calculated material properties 

Bl=LOG(SUl/SEl)/LOG(EUl*E/SEl) 'hardening exponent of PLATE mafl 
A1=SE1/((SE1/E) A B1) 'strength coeff. of mat’l 

B2=LOG(SU2/SE2)/LOG(EU2*E/SE2) 'hardening exponent of WELD mafl 
A2=SE2/((SE2/E) A B2) 'strength coeff. of mafl 

EK2=EU2/3 

SK1=A1*EK2 A B1 

SK2=A2*EK2 A B2 

BB1=L0G((SU1-SE2)/(SK1-SE2)) 

B B 1 =B B l/LOG((EU 1 -EE2)/(EK2-EE2)) 

AA 1=(SK1-SE2)/((EK2-EE2) A BB 1 ) 
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BB2=LOG((SU2-SE2)/(SK2-SE2)) 

BB2=BB2/LOG((EU2-EE2)/(EK2-EE2)) 

AA2=(SK2-SE2)/((EK2-EE2) A BB2) 

'poisson's ratio 

UXU1=.5-.2*SU1/E*((A1/SU1) A (1/B1)) 'poisson's ratio OF PLATE at ULT stress 
UXEL1=.3 

UXU2=.5-.2*SU2/E*((A2/SU2) A (1/B2)) 'poisson’s ratio OF WELD at ULT stress 
UXEL2=.3 

F 1 =LOG(UXU 1/UXEL 1 )/LOG((SU 1/SE1) A (1/B 1)) 

G 1 =. 3/((SE 1/A 1 ) A (F 1/B 1 )) 

F2=LOG(UXU2/UXEL2)/LOG((SU2/SE2) A (l/B2)) 

G2=.3/((SE2/A2) A (F2/B2» 

UXK=.3+(UXU2-UXEL2)/3 

SKU1=A1*(UXK/G1) A (B1/F1) 

SKU2=A2*(UXK/G2) A (B2/F2) 

FFl=LOG((UXU 1 -UXEL2)/(UXK-UXEL2)) 

FF1=FF1/L0G((SU 1 -SE2)/(SKU 1 -SE2)) 

GG 1 =(UXK-UXEL2)/((S KU 1 -SE2) A FF1) 

FF2=LOG ((UXU2 -UXEL2)/(UXK-UXEL2)) 

FF2=FF2/LOG((SU2-SE2)/(SKU2-SE2)) 

GG2=(UXK-UXEL2)/((SKU2-SE2) A FF2) 


R1S=0: R1F=1! : NR1=41 : DR1=(R1F-R1S)/(NR1-1) 'Y=0 to 1 

FOR 1=1 TO NR1 

R1=R1S+(I-1)*DR1 ' R!=Y/H 

YlsRl^Hl y/L=%h/L 

GOSUB T2SUB 

WRITE #24ll,T2(l),T2(2).T2(3),T2(4) 

PRINT R 1 ;T2( 1 );T2(2);T2(3);T2(4) 

NEXT I 

CLOS T #2 

STOP 

T2SUB: 

SX=0 

SY=0 

SXY=0 

SK=0 

FOR M=1 TO Ml 

A=M*3.1416 

AH2=A*H2 

AX2=AH2 

AH1=A*H1 

AY1=A*Y1 

EY1=EXP(AY1) 

EYN1=EXP(-AY1) 

SHY 1=(EY 1 -E YN 1 )* .5 

CHY1=(EY1+EYN1)*.5 

EH1=EXP(AH1) 

EH1N=EXP(-AH1) 
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SHH1=(EH1-EH1N)*.5 

CHH1=(EH1+EH1N)*.5 

E2H=EXP(2*AH1) 

E2HN=EXP(-2*AH1) 

SH2H1=(E2H-E2HN)*.5 

C1=SIN(AH2)*EH1N 

CX=(H1-1/A)*CHY1-Y1*SHY1 

CX=CX*C1*C0S(AX2) 

SX=SX+CX 

CY=Y1*SHY1-(H1+1/A)*CHY1 
CY =C Y*C 1 *COS (AX2) 

SY=SY+CY 

CXY=H1*SHY1-Y1*CHY1 
CXY=CXY*C 1 * SIN(AX2) 

SXY=SXY+CXY 

CK=H 1 *CHH 1 -(H 1 +2/ A)* S HH 1 
CK=CK*C 1 *COS( AX2)/AH 1 
SK=SK+CK 

NEXT M 

K2=1/(H2+4*SK) -EQUATION V-16 

Q 1= 1 -(GG 1 * ( AA2 A ( 1 /B B 2)))/(GG2* ( AA 1 A (1/BB1))) 

Q 1=Q 1 *K2*GG2/(AA2 A ( 1/B B2)) 

Q1=Q1/(1+(A2 A (1/B2))/(A1 A (1/B1))) 

K3=A2*Q1 A B2 -EQUATION V-19 

SSX=4*SX*K3 

SSY=(H2+4*SY)*K3 

SSXY=4*SXY*K3 

K4=(FF2*BB2+1)*B2/BB2 

K5=((SU2-SE2) A K4)/SU2 

T2(1)=SSX*K5 -equs V-22 

T2(2)=SSY*K5 

T2(3)=SSXY*K5 

T2(4)=(1+(T2(1)) A 2+T2(2) A 2-(1+T2(1))*T2(2)+3*T2(3) A 2) A .5 'equ V-23 

RETURN 

END 
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APPENDIX C 


MATERIALS ANALOG 


There are three noble disciplines from which mechanical behavior of most materials may be 
modeled: elasticity, viscosity, and plasticity. Their responses are unique to their content of force, dis- 
placement, and time variables. These three basic disciplines may be combined (fig. C-l ) to further define 
mechanical behavior of more complicated, though common, engineering materials. 

Viscoelasticity is emerging into a discipline of its own. It has borrowed from the large body of 
elasticity principles and coupled them with time-displacement behavior of fluids to represent polymer 
material behavior. 


Bingham materials are represented by visco-plastic behavior. Most common of these materials 
are wet paints and printer’s ink. Plastic-elastic behavior is inherent to most polycrystalline materials 
stressed beyond the elastic limit. 
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Figure C-l. Mechanical representation of common engineering materials. 
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